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In situ DNA barcoding of Trichoderma 
in soil reveals a narrow community of 
opportunistic species § 

Martina A. Friedl & Irina S. Druzhinina* 


Abstract: In this paper we report on the in situ diversity of the mycotrophic fungus Trichoderma (teleomorph 
Hypocrea, Ascomycota, Dikarya) revealed by a taxon-specific metagenomic approach and oligonucleotide 
DNA barcoding. We designed a set of genus-specific ITS1 and 2 rRNA primers and constructed a clone library 
containing 411 molecular operational taxonomic units (MOTUs). The overall species composition in soils of the 
two distinct ecosystems of the Danube floodplain consisted of 15 known species and 2 potentially new taxa. The 
latter taxa accounted only for 1.5 % of all MOTUs, suggesting that almost no hidden or uncultivable Hypocrea/ 
Trichoderma species is present at least in these temperate forest soils. The species were unevenly distributed 
in vertical soil profiles although no universal factors controlling the distribution of all of them (chemical soil 
properties, vegetation type, and affinity to rhizosphere) were revealed. Our data suggest that only a relatively 
small portion of Hypocrea/Trichoderma species is adapted to soil as a habitat and that the interaction between 
these species should be considered in a screening for Hypocrea/Trichoderma agent(s) of biological control of 
pests. 

Zusammenfassung: Die vorliegende Publikation zeigt Ergebnisse bezuglich der in-situ Diversitat des 
mycotrophen Pilzes Trichoderma (teleomorph Hypocrea, Ascomycota, Dikarya) basierend auf einer 
artspezifichen, metagenomischen Studie. Dazu wurden gattungsspezifische rRNA Primer fur ITS 1 und und ITS2 
entwickelt und damit eine umfassende Klonbibliothek (411 molecular operational taxonomic units = MOTUs) 
angelegt. Die Zusammensetzung der Arten in den beiden Bodenproben aus zwei verschiedenen Okosystemen 
der Donauauen umfasste 15 bekannte Arten und zwei potentiell neue, wobei letztere nur 1,5 % alle MOTUs 
ausmachten. Daraus lasst sich ableiten, dass zumindest in Waldboden der gemaBigten Zone keine versteckten 
Oder nicht kultivierbaren Arten von Hypocrea/Trichoderma existieren. Die verschiedenen Arten waren 
unterschiedlich in den vertikalen Bodenprofilen verteilt. Es konnte allerdings kein Zusammenhang zwischen 
globalen Einflussfaktoren wie den chemischen Eigenschaften des Bodens, der umliegenden Vegetation Oder der 
Affinitat zur Rhizosphare nachgewiesen werden. 
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INTRODUCTION 

Trichoderma species are one of the 
most frequently isolated conidial fungi 
(teleomorph Hypocrea, Hypocreales, 
Ascomycota, Dikarya). They are 
known from an innumerable diversity 
of natural and artificial substrata what 
demonstrates their high opportunistic 
potential and adaptability to various 
ecological conditions (Druzhinina et 
al. 2011). The exploitation of the lat¬ 
ter properties of Trichoderma in bio¬ 
technology and agriculture made the 
genus well studied and brought it in 
focus of numerous -omic studies [see 
Lorito et al. (2010) for the review] in¬ 
cluding the three complete genome 
sequences recently released for pub¬ 
lic access - H. jecorina (Martinez et al. 
2008), H. virens and H. atroviridis (Kubicek 
et al. 2011). The comparative analy¬ 
sis of these genomes with each other 
and other Ascomycota revealed that 
the outstanding antagonistic ability of 
Trichoderma spp. against plant patho¬ 
genic fungi (known as necrotrophic hy¬ 
perparasitism or mycoparasitism) is the 
ancestral state for the genus (Kubicek et 
al. 2011). Several studies have docu¬ 
mented that Trichoderma harzianum 
sensu lato, T. asperellum and T. asper- 
e/loides are highly rhizosphere com¬ 
petent and are capable to stimulate 
growth and immune defense of plants 
(reviewed by Harmann et al. 2004). 

The current diversity of the holomor- 
phic genus Hypocrea/Trichoderma 
is reflected in approximately 160 spe¬ 
cies (Samuels 2005; Druzhinina et al. 2006; 
Kubicek et al. 2008, NCBI Taxonomy 
Browser for February 2011), the majority 
of which have been recognized on the 
basis of DNA sequence analysis and 
molecular phylogeny of pure cultures 
and/or herbaria specimens. A multi¬ 
gene phylogeny established in several 
studies [see Kubicek et al. 2008 for re¬ 
view and Jaklitsch (20011) for updates] 
allowed the development of reliable 
tools for molecular species identifica¬ 
tion based on DNA barcoding, i.e. the 
analysis of DNA sequence polymor¬ 
phism in the internal transcribed spac¬ 


ers 1 and 2 (ITS1 and 2) of the rRNA 
operon and other loci (Druzhinina et al. 
2005; Kopchinskiy et al. 2005, www.lSTH. 
info). 

The complete inventory of te¬ 
leomorph forming species in Central 
Europe resulted in 75 species (Jaklitsch 
2009, 2001). In addition to it about 
a dozen of anamorphic species (for 
which no sexual stages are known) 
were reported in surveys of European 
soils and in taxonomic studies of 
Wuczskowski et al. 2003; Kraus et al. 2004; 
Samuels et al. 2006b; Jaklitsch et al. 2006; 
Komon-Zelazowskay et al. 2007; Hagn et 
al. 2007; Druzhinina et al. 2008; Migheu et 
al. 2009; Zachow et al. 2009 and Meincke 
et al. 2009. Thus, the total number of 
species Hypocrea/Trichoderma so far 
known in Europe may be about one 
hundred. 

The major ecological niche for 
Hypocrea/Trichoderma may be de¬ 
duced from the distribution of holomor- 
phic species (when both teleomorphic 
and anamorphic stages observed), 
which constitute the major genetic 
pool for the genus: they are associated 
with dead wood on different stages of 
its decay and with sporocarps of oth¬ 
er fungi (Jaklitsch 2009). Trichoderma 
species were initially believed to be 
among the dominant taxa inhabit¬ 
ing soil ecosystems (Klein and Eveleigh 
1998). However this statement was not 
supported by the application of high- 
throughput sequencing methods to 
study the cultivation-independent fun¬ 
gal diversity in soil. These studies did not 
reveal any abundance of Trichoderma 
listing it among minor taxa (Buee et al. 
2009; Lim et al. 2010). 

Thus, molecular ecology and ge¬ 
nomics of the genus indicate that pres¬ 
ence of Hypocrea/Trichoderma in soil 
where it may be either a saprotroph 
or establish various associations with 
plants and animals (biotrophy) may 
be driven by the general mycotrophy 
including various forms of mycopara¬ 
sitism, combined with broad environ¬ 
mental opportunism (Rossman et al. 
2003; see Druzhinina et al. 2011 for the 
review). 


The aim of this research was to 
use a taxon-specific metagenomic 
approach and oligonucleotide DNA 
barcode available for Hypocrea/ 
Trichoderma to explore the diversity of 
the genus in a soil profile and to identify 
the factors which control the size/oc¬ 
currence of its infrageneric communi¬ 
ties. 


METHODS 

Sampling sites 

Two sampling sites were chosen 
in the River Danube National Park 
“Nationalpark Donau-Auen” (Austria). 
They were located a few hundred me¬ 
ters apart representing the essentially 
different biotopes (Table 1). The beech 
forest site (+48° 9’ 28 N, +16°32'19 W, 
alt. 162 m) belongs to the hard wood 
riparian forest, which is situated above 
the seasonal Danube overflow level 
(approximately 3 m above the water 
level). The aspen forest site (+45° 30’ 43 
N, +73° 32’ 44 W, alt. 141 m) represents 
the softwood riparian forest regularly 
flooded during the seasonal Danube 
overflows. Detailed botanical descrip¬ 
tion of both ecosystems is given in 
Table 1. 


Soil sampling and detec¬ 
tion of soil properties 

In both ecosystems, pits of about 
100 cm of depth were made and soil 
samples of about 200 g were taken 
from each soil horizon using a sterile 
knife. The soil samples were immediate¬ 
ly separated from roots and large par¬ 
ticles, air dried and sieved through the 
metal mesh to collect fine particles < 
2 mm (“fine earth fraction”). Thereafter 
each soil was spread into a sterile tray 
and divided into 2x2 equally sized 
fragments, two of which were discard¬ 
ed while the remaining two were thor¬ 
oughly mixed and again spread on the 
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Table 1. Ecosystem description and soil properties. 


Beech forest site 


Aspen forest site 


Types of ecosystem 
Flooding frequency 


Hard wood riparian forest patch 
rare to none 


Clearing of the soft wood riparian 
forest 

regular 


Botanical description 


Trees 

dominant 

rare 

Shrub and regeneration 

dominant 

rare 

Herbs 

dominant 

rare 

Fallen trees 


tall, thick 

Ulmus laevis, Populus canescens, Fraxinus 
excelsior 

P. nigra, P. alba, Salix rubens 
lighted 

Crataegus monogyna 
Populus sp., Acer sp., Salix sp. 
lighted 

Aegopodium podagraria 
Salvia glutinosa, Poaceae, Cyperaceae 
abundant 


moderately tall, thin 
P. alba, P. nigra, P. canescens 
S. rubens 
dense 

Rubus idaeus 

Quercus robur, Populus sp. 
dense 

Allium ursinum, Ae. podagraria 

Carex acuta, Urtica dioica, 
Epilobium sp.. Salvia glutinosa, 
Solidago gigantea 

none 


Soil properties 


Litter layer 



thick, well developed 

thin, undeveloped 

Soil type (FAO) 



cambic fluvisol 


calcaic fluvisol with high level of 
sandy sediments 

Soil moisture content 


moderate 



high 


Siltation 



moderate 



high 


Vertical profile 



well differentiated 


well differentiated 


horizon name 

AO 

A 

B 

C 

A 

BC 

C 


thickness (cm 
of 1 m section) 

5 

10 

40-50 

20+ 

10 

20-30 

50-60 


soil color 

very 

dark 

brown 

dark 

grayish 

grayish brown 

very dark 
grayish 
brown 

light olive brown 


soil color code 
(Munsel) 

1OYR 2/2 

10YR 

4/2 

2.5Y 5/2 

2.5Y3/2 

2.5Y 5/3 


relative root 
density 

no roots 

very 

high 

high 

moderate 

high 

low 

no roots 

Chemical 

characteristics 

pH 

7.77 

7.90 

8.10 

8.24 

7.60 

8.15 

8.23 


C,% 

8.16 

2.38 

1.82 

2.50 

6.52 

1.04 

2.51 


N, % 

0.37 

0.15 

0.06 

0.04 

0.43 

0.07 

0.04 


C/N 

22.05 

15.87 

30.33 

62.50 

15.16 

14.86 

62.75 


The table was originally published in Friedl & Druzhinina 2012, the special issue of Microbiology “Trichoderma - from 
Basic Biology to Biotechnology”, January 2012 


STAPFIA 96 (2012) | 181 


































Friedl & Druzhinina • DNA barcoding of Trichoderma STAPFIA 96 (2012): 179-193 


Table 2. Genus and species-specific PCR primers designed in this study. 

Primer 

Locus 

Application 

Suggested 

pair 

Specificity 

Direction 

Primer sequence 

Trirevl 

ITS2 

ciPCR 

ITS5 

genus, Hypocrea 
and Trichoderma 

reverse 

5'-C ATTTC (A/C) GAAAGTTGGGGTG-3' 

Trirev2 

ITS2 


ITS5 



5'-CATTTC(A/C)GAAGTTTGGGGTG-3' 

Trirev3 

ITS2 


ITS5 



5'-C ATTTC (A/C) GAAAGTTTGGGT G-3' 

Trirev4 

ITS2 


ITS5 



5'-CATTTC (A/C) GAAAGTTGGGTG-3' 

Trirev5 

ITS2 


ITS5 



5'-CATTTC (A/C) GAAGTTGGGTG-3' 

Trirev6 

ITS2 


ITS5 



5'-C ATTTC (A/C) GAAGTTTGGTG-3’ 

citrol_64 

tefi 

ciPCR and qPCR 

LLErev 

species, 

H. schweinitzii 

forward 

5'-CGCTACTGCCTTCAGACCAC-3' 

asp7_60 

tefl 


LLErev 

species, 

T. asperellum 

forward 

5'-GCTTTGCCAGTCTACCTACC-3' 


The table was originally published in Friedl & Druzhinina 2012, the special issue of Microbiology “Trichoderma - from 
Basic Biology to Biotechnology", January 2012 


same tray for subsequent subsampling 
(Robertson 1999). Finally 50 g of each soil 
was stored at -20° C for further molecu¬ 
lar and chemical investigations. 

Soil classification and the detection 
of soil horizons were done directly at 
the two sampling spots. The soil color 
was defined using a standard color 
scale for soil science (Munsell Soil Color 
Charts, U.S. Dept, of Agriculture). The 
horizons were visually distinguished 
based on their color and later on dif¬ 
ferentiated based on their chemical 
properties. All chemical analyses were 
performed using the fine earth fraction. 
To measure pH, 1 g of soil was suspend¬ 
ed in 100 ml of 1M KCI and shaked for 
one hour. Finally, pH was determined 
with a glass electrode. The total nitro¬ 
gen content was determined accord¬ 
ing to the Kjeldahl method (see Batjes 
1996 for a reference) on a Vapodest 30 
(Gerhardt, Germany). The total organ¬ 
ic carbon content was measured using 
the Liechtenfelder method (see Batjes 
1996 for a reference), which oxidized 
carbon with potassium dichromate, 
and quantifies the generated Cr 3+ pho¬ 
tometrically (DIN 19684). 


Development of genus- 
specific primers 

The genus-specific primers were 
developed based on the master align¬ 
ment of ITS 1 and 2 from 88 reference 
strains of Hypocrea and Trichoderma 
(Druzhinina et al. 2005) complemented 
by the new species described since 
that time (Jaklitsch et al. 2005; Jaklitsch 
et al. 2006a; Jaklitsch et al. 2006b; 
Overton et al. 2006; Samuels et al. 2006b; 
Komon-Zelazowska et al. 2007; Jaklitsch et 
al. 2008; Jaklitsch 2009, 2011). All prim¬ 
ers are reverse primers and are com¬ 
plementary to the forward primer ITS 5 
(5‘-GGA AGT AAA AGT CGT AAC AAG 
G-3‘) (Table 2). The position of primers is 
shown on www.ISTH.info/metagenom- 
ics. The verification of primer specific¬ 
ity was done with reference cultures 
presented in Druzhinina et al. (2005). 
The general test of possible unspecific 
annealing was done by sequence 
similarity search against NCBI GenBank 
(May 2006) adjusted for short query se¬ 
quences. 


Isolation of DNA and 
ciPCR 

A representative subsample of 1 g 
of fine soil fraction was thawed and 
dried over night in a drying chamber 
at 90°C. Afterwards, the individual soil 
samples were well homogenized and 
DNA was extracted using a FastDNA 
Spin kit for soil (MP Biobiomedicals, 
Germany). The DNA was then used as 
a template for ciPCR (culture indepen¬ 
dent PCR) with Hypocrea/Trichoderma 
specific rRNA primers as listed in Table 
2. CiPCR products of all six primers 
were combined and purified using 
QiaQuick PCR purification kit (Qiagen, 
Germany). 

CiPCR reactions were carried out 
in a total volume of 50 pi containing 
2.5 mM MgCI 2 , 10 mM Tris- HCI pH 9.0, 
50 mM KCI, 0.1 % (v/v) Triton X-100, 0.4 
pM of each primer, 0.2 mM of each 
dNTP and 0.5 units of Taq-Polymerase 
(Promega, Madison, Wl). The ampli¬ 
fication program consisted of: 1 min 
initial denaturation (94°C), 30 cycles of 
amplification (1 min 94°C, 1 min 52° C, 
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IT SI and 2 of the rRNA gene cluster as 
the universal DNA barcode markerfor 
fungi: 

• multiple copiesin the genome 

• easyto amplifyfrom pure cultures and 
envrronm ental sam pies 

• containsconservedandhypervaiiable 
areas 

■ resolves majority of fung atsubgeneric 
levels (clades or species) 


Requirements of TrickOKty (DRUZHININA et al. 2005), an online tool for oligonucleotide DNA 
barcode of the frequent species of Hypocrea/Trichoderma; 

• m olecular evolutron and diversity of the group hasbeen studied 

• auniversalbarcodemarkerhasbeen defined 

• genus-specific hallm arkshave been established 

• a database of species-specific (a diagnostic combination ofup to J ohgpnucleotide DNA barcodes 
per each species) hasbeen establishedbasedon known infraspecific diversity of the most comm on 
species 

Advantage an absolute identification without custom evaluation of sequence similarity values 
Limitations available for the most frequent species with known infraspecific variability 


ITS1 

diagnostic area 


5.8 S rRNA gene 
none-diagnostic area 


ITS2 

diagnostic area 


mi mill ilium i ini iiiiiiiiiiiiiiiii rii niimiiii ihiiii ini m iiiiiiiiiiiiii m 


MM 

1_ 


tan 


ina 

i 


During sequence processing the query DNA sequence is screened for the presence of five genus-specific oligonucleotide barcodesthatshouldbe 
arrangedin agven order. If positive, the sequence is attnbutedtothe target group of organism s (genus Hypocrea/Trichoderma) After that the 3 
regions for the search of species-specific DNA barcodes are assisted 


Region I ^_ Region2 

an iiiniaiufliitainiiaiii 

1 _ 


in in mi in 
■a 


in hi mil 


Region3 


aiiai naim It 1 ’ 111 III III I 
I 1! 

_ I 


Identification isbased on the meaningfulcombination oflaxon-specific oligonucleotide DNA Barcode found. 


Fig. 1: Principle of the oligonucleotide DNA barcoding implemented in TrichO Key (Druzhinina et al., 2005) 


1 min 72°C) and a final extension pe¬ 
riod of 7 min at 72°C. 

Amplification products were visual¬ 
ized on 1% TAE agarose gels and the 
ciPCR amplicons at a size of 600 kb 
were excised, repurified using QiaQuick 
Gel excision kit (Qiagen, Germany) 
and suspended in 20 pi water. 


Construction of clone 
libraries and sequencing 

The ciPCR products were subcloned 
using the pGEMTeasy (Promega) stan¬ 
dard procedure. At least 100 colonies 
were selected from each soil sample. 
Plasmids were extracted using a stan¬ 
dard miniprep method. In order to con¬ 
trol if the plasmid DNA was carrying the 
insert, the plasmid DNA was digested 
with 5u of the restriction enzymes Not] 
and EcoRl enzymes (Fermentas, Glen 


Burnie, USA). For positive samples 0.5 
pg of the plasmid DNA was finally used 
for automated sequencing (Eurofins 
MWG, Ebersberg, Germany) from both 
directions. 

Representative alleles for all spe¬ 
cies from each horizon have been 
deposited in the NCBI GenBank da¬ 
tabase (accession numbers. T. harzia- 
num GQ981644, GQ981645, GU060097 
- GU060105; T. cerinum GU060106, 
GU060107; T.pleuroticola GU060108, 
GU060109, H. alni GU060110; H. virens 
GU060111; T. sp. MOTU 1A 64 GU060112 

GU060117; T. brevicompactum 
GU060118 - GU060120; H. pachybasi- 
odes GU060121; uncult. H. pachypat- 
lida GU060122 - GU060125: T. longi- 
brachiatum GU060126 - GU060128; H. 
schweinitzii GU060129, T. asperellum 
GU060130 - GU060137; T. sp. MOTU 2B 
48 GU060138, GU060139. The align¬ 
ment matrix is available upon request. 


Species identification by 
DNA barcode and diver¬ 
sity assessment 

All sequences were aligned in 
GeneDoc 2.6 (Nicholas and Nicholas 
1997) using the guidance for 
Hypocrea/Trichoderma ITS1 and 2 
alignment provided by Druzhinina et 
al. (2005), available on-line at http:// 
www.isth.info/tools/master.php. For 
species identification the complete 
set of sequences was submitted to 
the DNA oligonucleotide barcode 
program TrichO Key (http://www.isth. 
info/tools/molkey/index.php; Druzhinina 
et al. 2005). All sequences contained 
the set of five genus-specific hall¬ 
marks and were therefore attributed 
to Hypocrea/Trichoderma. Hence, 
tests for chimerical sequences were 
not required. All unusual ITS1 and 2 al¬ 
leles have been further analyzed by 
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sequence similarity search against the 
NCBI GenBank, TrichoBLAST (Kopchinskiy 
et al. 2005) and the sequence data¬ 
base of the TUCIM the Collection of 
Industrial Microorganisms of the Vienna 
University of Technology that currently 
contains more than 4000 Hypocrea/ 
Trichoderma strains with 5500 core 
nucleotide sequences including ITS1 
and 2. Based on individual mismatches 
found in the otherwise conserved areas 
of ITS1 and 2 (for example, genus-spe¬ 
cific hallmarks or 5.8S rRNA gene, see 
Druzhinina etal. 2005), 20% of sequences 
have been diagnosed to contain sin¬ 
gle sequencing errors. Four sequences 
contained polymorphic sites in the di¬ 
agnostic regions of both ITS1 and 2 and 
therefore have been diagnosed as po¬ 
tentially new alleles (T. sp. MOTU 1A 64 
for Section Longibrachiatum and T. sp. 
MOTU 2B 48 for Section Trichoderma ). 
The principle of the oligonucleotide 
DNA barcode is shown if Fig. 1. 


Development and verifi¬ 
cation of species-specific 
tefl primers for qPCR 

In order to design species-specif¬ 
ic qPCR primers, representative se¬ 
quences of the 4 lh large intron of tefl 
gene coding the elongation factor 1 
alpha for the whole genus Hypocrea/ 
Trichoderma were retrieved from the 
multilocus database of phylogenetic 
markers ( www.ish.info ) and automati¬ 
cally aligned in Clustal X (Thompson et 
al. 1997). The representative tefl se¬ 
quence of a target species was sub¬ 
mitted to NCBI sequence similarity 
search tool (blastn) and all homolo¬ 
gous vouchered sequences (N>20) at¬ 
tributed to the same species retrieved 
and added to the initial alignment. 
Subsequently, species-specific diag¬ 
nostic regions were manually selected 
for each target species in a way that 
they contain the minimal level of infra¬ 
specific polymorphism. In parallel, the 
same alignments were used to design 
degenerate species-specific primers 
in HYDEN software (Linhart and Shamir 
2005). The annealing temperature and 


secondary structure of oligonucle¬ 
otides designed based on both ap¬ 
proaches have been estimated using 
Gene Runner (Gene Runner 3.0 soft¬ 
ware) and SMS PCR Primer Stats tools 
http://www.bioinformatics.org/sms2/ 
pcr_primer_stats.html. Specificity of se¬ 
lected primers was first tested against 
NCBI GeneBank database (automati¬ 
cally optimized for short queries) and 
then verified by PCR reactions (see 
Jaklitsch et al. 2006a for conditions) with 
reference DNAs from pure cultures of all 
genetically close and several members 
of the neighbor clades. The resulting 
potentially species-specific oligonucle¬ 
otides were tested for the annealing 
efficiency by applying serial dilutions 
of a target DNA extract. Selectivity of 
designed primers was verified by sub¬ 
cloning the PCR product obtained 
from soil DNA extract. Around 40 oli¬ 
gonucleotides have been synthesized 
to hit the diversity of the most frequent 
temperate Hypocrea/Trichoderma 
species but only two primers targeting 
T. asperelium and H. schweinitzii (Table 
2) showed high specificity, selectivity 
and appropriate efficiency (>80%) to 
be applied to environmental samples. 


Semi-quantitative PCR 
assessment 

Quantitative PCR amplification was 
carried out with the iQ 5 Real-Time PCR 
detection system (Bio-Rad) in a 25 |il re¬ 
action containing 12.5 |il iQ SYBR Green 
Supermix (Bio-Rad), each primer at a 
concentration of 250 nM and sample 
corresponding to an initial concentra¬ 
tion of 0.5 |il of total DNA. Amplification 
was carried out with the following PCR 
program: initial denaturation for 3 min 
at 95°C, followed by 45 cycles consist¬ 
ing of 95°C for 15 s, 54.0°C (tefl), re¬ 
spectively, for 20 s, and 72°C for 20 s. 
Successful amplification was verified by 
determination of the melting tempera¬ 
ture and by agarose gel electrophore¬ 
sis. For each species a series of dilutions 
was performed to assess the efficiency 
of the PCR. The results of the real-time 
PCR were analyzed with the iQ 5 opti¬ 
cal system software (Bio-Rad). Using the 


PCR base line subtracted mode, the 
threshold cycle was calculated for all 
samples and the amplification efficien¬ 
cy for each primer was determined. 


Statistical analysis 

Cultivation-independent diversity 
and species richness of Hypocrea/ 
Trichoderma was assessed in Estimates 
8.2. Rarefaction curves were comput¬ 
ed for each soil profile. The number of 
species was quantified for 100 random 
combinations of 1 to N sequences and 
also by performing 100 bootstrap pseu- 
doreplicates implemented in Estimates 
(Colwell 2005). 

Statistical analyses of metage- 
nomic data and in vitro infrageneric 
interactions were done in Statistica 
6.1 (StatSoft Inc., Tulsa, USA) software 
package using basic data exploration 
tools, correlation matrices, variance 
analyses and multifactorial techniques 
(factor analysis and cluster analysis). 
Discrete color plots representing the re¬ 
sults of pair wise interactions between 
different Hypocrea/Trichoderma strains 
were constructed on the basis of a re¬ 
sult of a two-way joining cluster analysis 
implemented in Statistica 6.1 with con¬ 
sequent reordering of both variables 
and cases in order to reflect phyloge¬ 
netic groups inside the sample. 


RESULTS 

The pedogenesis of both 
sites is strongly influenced 
by the river 

The two sampling sites representing 
essentially different biotopes (Table 1) 
in the riparian forest were selected in 
the River Danube National Park south¬ 
east of Vienna, Austria. We performed 
a 1 m deep vertical soil cut into both 
ecosystems in a place which is equally 
distanced from surrounding big trees 
and shrubs and which therefore does 
not represent an exclusive rhizosphere 
of any plant species. 
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Fig. 2: In situ diversity of Hypocrea/Trichoderma in vertical profiles of two soil cuts in a hard wood beech and a soft wood aspen forest 
sides respectively. Numbers in colored cycles indicate the number of MOTUs of each species recovered: red, yellow and blue cycles 
correspond to dominant, sub-dominant and rare species, respectively. The middle right insert shows the schematic profile of the flood 
plain ecosystem of the river Danube and images of the sampling sites. * - T. longibrachiatum - H. orientalis species pair. The figure was 
originally published in Friedl & Druzhinina 2012, the special issue of Microbiology "Trichoderma - from Basic Biology to Biotechnology", 
January 2012 http://mic.sgmjournals.Org/content/158/l/69.long. 


STAPFIA 96 (2012) | 185 























Friedl & Druzhinina • DNA barcoding of Trichoderma 


STAPFIA 96 (2012): 179-193 


20 


18 - 


16 - 
IA 

O 

a 14 - 

Q. 

i/i 

H- 

2 12 - 

a 

_Q 


c 


8 - 


6 - 


4 ■ i-1-1-1-1-1-1-1-1 

0 50 100 150 200 250 300 350 400 450 


computed number of individuals 

Fig. 3 : Rarefaction curve estimating the total Hypocrea/Trichoderma community diver¬ 
sity as a function of clone library size. The curve was computed from 100 bootstrap 
replicates (blue) and the expected richness function (red) is the number of species 
estimated (random sampling without replacement) from our sequences dataset. 
Vertical bars show the standard deviation. The figure was originally published in Friedl & 
Druzhinina 2012, the special issue of Microbiology “Trichoderma - from Basic Biology to 
Biotechnology”, January 2012 http://mic.sgmjournals.Org/content/158/l/69.long. 



The major difference between the 
two soil cuts was in the presence of the 
thick (min 5 cm) soil litter layer (horizon 
A0, FAO) in the beech site which was 
absent in the aspen forest. Instead, the 
latter site was covered by the dense 
herbal layer dominated by the blos¬ 
soming Alnus urcinum with several co¬ 
dominated herbal species (see Table 
1). Correspondingly the beech forest 
soil profile was evenly penetrated by 
roots from different plants while in the 
aspen forest a high root density was 
observed only in the upper horizon. 

The pedogenesis of both sites is 
strongly influenced by the river and re¬ 
sults in formation of fluvisols (FAO classi¬ 


fication: J) with a high level of calcaric 
soil material and sedimentation. The 
flooded soil in the aspen forest contains 
significantly more moisture compared 
to the beech forest (Table 1). Both sites 
were characterized by well developed 
soil profiles and presence of at least 
three horizons in FAO nomenclature: 
the soil of the beech forest was char¬ 
acterized by three clear horizons A, B 
and C; at the aspen forest the B layer 
contained little alteration products 
and was mixed with the parent mate¬ 
rial, therefore the horizons at this sam¬ 
pling spot were classified as A, BC and 
C (Table 1). The chemical properties of 
the soil are detailed in Table 1. 


DNA barcode revealed 
restricted diversity of 
Hypocrea/Trichoderma in 
temperate soil 

We designed the set of genus- 
specific ciPCR (culture independent 
PCR) reverse primers Trirev (1 to 6) 
5'-C ATTTCMG [A 2 /A 3 ] G [T 2 /T 3 ] [G 2 /G 3 /G 4 ] 
TG-3' (Tm 61 °C) (Table 2) in the 3 prime 
end of the internal transcribed spacer 
2 (ITS2) of the rRNA gene cluster (Table 
2). The position of primers is illustrated 
on http://www.isth.info/methods/ 
method. php?method_id=l 2. When 
applied together with ITS5 forward 
primer (White et al. 1991), these primers 
amplified a 510 - 540 bp fragment cov¬ 
ering the complete diagnostic area of 
ITS1 and 2 including all five genus-spe¬ 
cific hallmarks (Druzhinina et al. 2005). 

Application of ITS5/Trirev primer pair 
to ciPCR with DNA extracts from soils 
samples showed different results for in¬ 
dividual primers and also unequal effi¬ 
ciency of Trirev primers in different ho¬ 
rizons. The highest affinity of Trirev prim¬ 
ers was detected for the horizon A fol¬ 
lowed by the two other soil horizons of 
the beech forest site. PCR efficiency in 
soil of aspen forest was essentially lower. 
The pooled ciPCR products of ITS5/Trirev 
obtained for each soil horizon were 
used for subsequent clone libraries. 

In total 411 ITS1 and 2 rRNA mo¬ 
lecular taxonomic units (MOTUs) were 
recovered (Fig. 2). All sequences had 
diagnostic genus-specific hallmarks 
(Druzhinina et al. 2005) suggesting high 
selectivity of Trirev primers. All MOTUs 
were identified as 15 known species 
and two putatively new taxa (Table 3). 
The rarefaction curve suggested that 
the species richness was close to satu¬ 
ration as additional 100 MOTUs could 
reveal only one or two additional spe¬ 
cies (Fig. 3). 

MOTUs frequencies were used to 
calculate the diversity index as an in¬ 
direct quantitative measure of com¬ 
munity composition. The dominant 
species in both sites (54 and 28.6 % for 
beech and aspen forests respectively) 
was T. asperellum sensu stricto (Samuels 
et al. 2009). 
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Beech forest Aspen forest 




ct' 1 ct' 1 

Fig. 4: Results of the semi-quantitative PCR analysis based on species-specific tefl primers de¬ 
veloped for T. asperellum and H. schweinitzii respectively. As the efficiency of primers binding 
the environmental (soil) DNA was slightly different between these species only semi-quantita- 
tive comparison is possible, ct' 1 (the number of threshold cycles/1) corresponds to abundance 
of the target DNA in the sample. The figure was originally published in Friedl & Druzhinina 2012, 
the special issue of Microbiology “Trichoderma - from Basic Biology to Biotechnology", January 
2012 http://mic.sgmjournals.org/content/158/1 /69.long. 


The beech forest soil was co-dom- 
inated by T. harzianum sensu stricto 
(Druzhinina et al. 2010b), T. rossicum, H. 
'pseudoharzianum' (Druzhinina et al. 
2010b, widely known as T. harzianum 
complex) and H. pachypallida (Jaklitsch 
2011) which contributed 11,9, 8 and 6 
% of the total diversity at this location. 
Sequences of T. longibrachiatum - H. 
orientalis species pair (Druzhinina et al. 
2008) added 5 per cents and could be 
considered as a subdominant taxon 
(Fig. 2). T. pleuroficola (2.5%), H. sch¬ 
weinitzii (1.2%), T. brevicompactum 
(0.8%), H. pachybasioides (0.8%) and 
taxonomic units of a putative new spe¬ 
cies from Section Longibrachiatum T. 
ct. sp. nov. MOTU 1A 64 (0.8%) were 
rare. MOTUs of H. atroviridis and H. Vi¬ 
rens were recovered with a frequency 
of less than 0.5 %. 

A qualitatively similar but quanti¬ 
tatively different species composition 
was observed in the aspen forest (Fig. 


2): the dominant T. asperellum was fol¬ 
lowed by co-dominant H. pachypal¬ 
lida (23%), H. ‘pseudoharzianum’ (16%) 
and T. rossicum (15.5%). The subdomi¬ 
nant species H. alni (5%) and T. cerinum 
(3%) were significantly less present. T. 
harzianum s. s. (2.4%), T. longibrachia¬ 
tum - H. orientalis species pair (1.8%), T. 
pleuroficola (1.2%), T. brevicompactum 
(1.2%), T. sp. C.P.K. 2974 (I. Druzhinina, un¬ 
published data) (1.2%) and a putative 
new species from Section Trichoderma 
T. sp. nov. MOTU 2B 48 were rare; the ITS 
sequence of H. schweinitzii was recov¬ 
ered only once. 

Thus, among 17 recovered taxa, we 
detected nine species present in both 
sites and four species unique for each 
ecosystem. In total, soils in both forest sites 
supported the community of at least 13 
co-existing Hypocrea/Trichoderma spe¬ 
cies although these communities were 
unequal in two locations. 


Highest Hypocrea/ 
Trichoderma diversity was 
detected in the moist ho¬ 
rizon of aspen forest soil 
that has a low root density 
and poor carbon content 

The analysis of Hypocrea/Tricho¬ 
derma species distribution in the verti¬ 
cal soil profile showed that T. asperel¬ 
lum was the most abundant in the soil 
litter layer of the beech forest (82 %) 
and also dominated in A and B hori¬ 
zons (53 and 59 % respectively) while in 
the horizon C it was only the third most 
abundant taxon (16 %). A similar distri¬ 
bution pattern was observed in the as¬ 
pen forest soil (Table 3). The correlation 
analysis showed a significant nega¬ 
tive correlation of T. asperellum abun¬ 
dance with pH values (Table 3) that 
decreases with depth, indicating that 
this species is associated either with the 
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soil litter layer or with the upper organic 
horizon A. 

The distribution of the next most fre¬ 
quent species H. pachypallida in the 
vertical soil profile is reverse proportion¬ 
al to that of T. asperellum (Fig. 2, Table 
3): in the beech forest soil it dominated 
the deepest C horizon (30 %) and was 
not detected in others while in the as¬ 
pen forest it was found in all horizons 
but with increasing frequency from top 
to down. The significant correlation of 
its occurrence with the horizon (Table 
3) but not with other factors (carbon 
content, pH and depth) suggests the 
presence of yet another parameter 
which controls the distribution of this 
species and which was not monitored 
in this study. 

T. rossicum - one of the co-dominant 
species in both locations - is also more 
frequent in the deeper soil horizons 
showing a significant positive correla¬ 
tion with the C/N ratio (Table 3). MOTUs 
of the T. longibrachiatum - H. orientalis 
species pair were the most frequent in 
the A horizon of the beech soil profile 
which has the highest density of roots 
(r= 0.8, P<0.05). 

No correlations were detected for 
H. ‘pseudoharzianum’ and T. harzia- 
num s. s. 

The highest diversity of Hypocrea/ 
Trichoderma (12 coexisting species) 
was found in the moist BC horizon 
of aspen forest soil, which has a low 
root density and the lowest carbon 
content of this study. The next most 
diverse infrageneric community (10 
species) was in the topsoil horizon A in 
the beech forest, which has the high¬ 
est root density. As these communities 
share six species but have no similar 
frequencies, consequently, no sig¬ 
nificant correlations between species 
richness and soil properties were de¬ 
tected. 

Since the genetic diversity of 
Hypocrea/Trichoderma in individual 
soil profiles appeared unexpectedly 
high we applied several statistical 
techniques to reveal possible cor¬ 
relations between the distributions 
of MOTUs from different species. The 
product-moment correlation analy¬ 
sis confirmed the significant (P<0.05) 


negative correlation between T. 
asperellum on one hand and H. 
pachypallida and T. rossicum on the 
other. No other significant correla¬ 
tions between dominant (co- and 
sub-) species were revealed suggest¬ 
ing that these species have different 
responses to the microecological 
conditions. 


Quantitative PCR confirms 
the vertical distribution of 
7. asperellum 

In order to have a second means 
of testing species distribution in soil pro¬ 
file we designed species-specific prim¬ 
ers for quantitative PCR (qPCR) based 
on the polymorphic fragment of the 
tefl gene. Although 40 primers were 
tested (see Methods for details), only 
those for T. asperellum and H. schwein- 
itzii showed to be highly selective and 
efficient when applied to soil DNA ex¬ 
tract (Table 2, Fig. 4). Using them, we 
detected a higher abundance of T. as¬ 
perellum and a relatively small amount 
of H. schweinitzii DNA in all seven 
samples, which is in accordance with 
the metagenomic data. It was also 
possible to detect a reduction of T. as¬ 
perellum biomass with depth while the 
distribution of H. schweinitzii remained 
constant. 


DISCUSSION 

To study the generic community 
of Hypocrea/Trichoderma in soil we 
selected a riparian forest in the River 
Danube National Park ("Nationalpark 
Donau-Auen”) southeast of Vienna, 
Austria. This ecosystem is unique, as it 
still resembles the original European 
river floodplain landscape free from 
anthropogenic loadings (Wuczskowski 
et al. 2003). The forest is characterized 
by high biodiversity including several 
otherwise rare plant species (Tockner et 
al. 1998). 


Temperate soil supports a 
limited diversity of highly 
opportunistic Hypocrea/ 
Trichoderma spp. 

Although the biodiversity of Higher 
Fungi is considered to be largely un¬ 
known (Hawksworth 1991), and studies 
using cultivation-independent meth¬ 
ods should thus result in the identifi¬ 
cation of a high percentage of un¬ 
known taxa, this is not the case with 
Hypocrea/Trichoderma. Our data sug¬ 
gest that at least in soils of temperate 
climate there is almost no hidden diver¬ 
sity of Hypocrea/Trichoderma: among 
411 MOTUs 407 were safely attributed 
to 15 existing species or to putatively 
new taxa that have previously been 
sampled. The diversity of Hypocrea/ 
Trichoderma in Europe consists of at 
least 75 holomorphic species (Jaklitsch 
2009) and 10-20 anamorphic species 
(see refs, in Introduction), in summary 
approaching one hundred taxa. Our 
finding of only a minor portion of po¬ 
tentially expected diversity (roughly 
15 %) is in agreement with the previ¬ 
ous hypothesis that soil itself is not the 
primary ecological niche for the ge¬ 
nus (Druzhinina et al. 2011). The similar 
outcome was also obtained by the 
pioneering metagenomic studies of 
HypocreaTtrichoderma: Hagn et al. 
(2007) used ITS1 fragment and found 
only about a dozen species in arable 
soil. Meincke et al. (2010) used a par¬ 
tial ITS1 and 2 and detected about 
20 Hypocrea/Trichoderma taxa in rhi- 
zosphere of Solanum tuberosum, al¬ 
though no undoubtful species identi¬ 
fication was made. This view is also in 
agreement with studies that applied 
high-throughput sequencing to reveal 
the actual in situ diversity of high fun¬ 
gi in soil: in these studies, Hypocrea/ 
Trichoderma MOTUs were found only 
at minor portions compared to other 
dominating groups of Ascomycota 
(Buee et al. 2009; Lim et al. 2010). 

The perception of Hypocrea/ 
Trichoderma as a common soil fungus 
is based on the abundant isolations 
from soil samples world-wide. However 
the qualitative analysis of the diver¬ 
sity revealed in such samples shows 
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the dominance by the same 15-20 
highly opportunistic species such as 
T. asperellum, T. aspereltoides, T. cf. 
harzianum, T. hamatum, T. atroviride, T. 
virens, T. longibrachiatum, T. gamsii, T. 
citrinoviride, T. koningiopsis, T. spirale, T. 
koningii complex etc. (see http://www. 
isth.info/materials/topic.php?material_ 
id=42 for details) which likely obtained 
the ability to saprotrophic growth in 
soil due to their general outstanding 
opportunistic potential as suggested 
based on genomes of T. atroviride and 
T. virens (Druzhinina et al. 2011; Kubicek et 
al. 2011). Moreover the antifungal ac¬ 
tivity of Hypocrea/Trichoderma spp. 
favors their detection, as they are able 
to suppress other fungi. The view of soil 
as the main determinant of Hypocrea/ 
Trichoderma ecological niche does 
not find its support. 

The qualitative composition of 
Hypocrea/Trichoderma community 
reveals that soil is inhabited by highly 
opportunistic species with cosmopoli¬ 
tan distribution as all taxa, except the 
3 MOTUs of putatively new species, are 
common and known from multiple iso¬ 
lates from numerous substrata (includ¬ 
ing soil) from temperate ecosystems 
world-wide. 


The occurrence of 
Hypocrea/Trichoderma 
species in different soil 
profiles is not determined 
by root density 

The PCR efficiency with Hypocrea/ 
Trichoderma- specific primers was es¬ 
sentially higher in the moist beech for¬ 
est (all horizons except the litter layer 
AO) compared to aspen site, while 
there was no visual difference between 
both sites when general fungal prim¬ 
ers were applied (data not shown). 
This finding indirectly suggests that the 
beech forest contains more Hypocrea/ 
Trichoderma biomass compared to the 
aspen soil. 

Moreover, not all species known 
to be abundant in this region have 
been detected in soil profiles. For in¬ 


stance, the most frequent teleomor- 
phic Hypocrea/Trichoderma species 
in Central Europe - H. minutispora 
(Jakutsch 2009) - was not found in our 
study at all. Despite of this, H. minutis¬ 
pora was abundantly found in air sam¬ 
ples in nearly the same region (within 
a few km 2 ) by the same methodolo¬ 
gy as here (M. A. Friedl, I.S. Druzhinina, 
unpublished), thus supporting the as¬ 
sumption that there is no methodical 
bias in our experimental procedures. 
Similar, many other very common lo¬ 
cal Hypocrea species (H. viridescens, 
H. rufa, H. puivinata, H. strictipilosa 
etc., Jakutsch et al. 2006; Overton et al. 
2006) were not detected in Danube 
floodplain soils indicating other eco¬ 
logical niches for their anamorphic 
stages. 

The observed distribution of MOTUs 
in the vertical soil profile also suggests 
the existence of ecological factors 
that determine the proliferation of 
Trichoderma. T. asperellum, which was 
the dominant species in this study, was 
reproducibly associated with upper 
soil horizons, but its abundance did not 
correlate with the carbon and nitrogen 
content of the soil or pH of soil solution. 
Also, the highest number of MOTUs for 
T. asperellum was found in the root-free 
A0 litter layer of the beech site, and we 
therefore suspect that the rhizosphere 
is not its prime habitat. However, it 
might be that the species is following 
other fungi highly abundant in soil litter 
and upper soil horizons. In contrast to 
that, H. pachypallida and T. rossicum 
were almost exclusively detected in 
the deepest soil horizons character¬ 
ized by the sufficient organic carbon 
and the nearly complete absence of 
plant roots. It suggests predominantly 
saprotrophic nutrition for the later two 
species in deep soil horizons. Moreover, 
these species should be capable of 
growing well under conditions of ni¬ 
trogen starvation as the amount of 
nitrogen in their habitat is one order 
of magnitude lower than that in the 
surface soil. Both the high species rich¬ 
ness of the deepest mineral soil hori¬ 
zons (8 species in each sampling site) 
and their unique species compositions 
compared to upper soil layers suggest 
competitive relations between tribal 


relatives (and/or with other myco- and 
microbionts) rather than their asso¬ 
ciations with abiotic characteristics of 
these soil horizons. 

The limited diversity of Hypocrea/ 
Trichoderma in studied soils may be 
attributed to the relatively high pH val¬ 
ues (around 8) which could potentially 
prevent the development of other 
species. However the same cosmo¬ 
politan and opportunistic species were 
also dominated mole acidic soils on 
Sardinia (pH =~ 5, Migheli et al. 2009) 
and were present in rhizospheric soil of 
Coffea arabica in Ethiopian highland 
forest (pH =~ 5.5, Mulau et al. 2010) sug¬ 
gesting that this is not the controller of 
the diversity. 


In situ diversity confirms 
the sympatric speciation 
within T. harzianum sensu 
lato and related taxa 

In most previous studies that used 
cultivation-dependent methods to 
quantify Hypocrea/Trichoderma in 
various habitats, T. harzianum sensu 
lato represented the most dominantly 
occurring species (Druzhinina et al. 2005; 
Migheli et al. 2009; Zachow et al. 2009 
and Druzhinina et al. 2010b for more ref¬ 
erences). Also in this study, we detect¬ 
ed a remarkable diversity of geneti¬ 
cally sibling species from Harzianum 
- Catoptron Clade (Chaveri et al. 2003; 
Druzhinina et al. 2010b) in nearly all soil 
samples. Interestingly, in all mineral soil 
horizons of the beech forest site T. har¬ 
zianum s.s. coexisted with a member of 
H. ‘pseudoharzianum' (Druzhinina et al. 
2010b) while in the aspen forest site the 
later one was found together with other 
members of the Harzianum - Catoptron 
clade: T. cerinum, T. pleuroticola and 
H. alni. This finding suggests the sym¬ 
patric speciation of T. harzianum s. s. 
and H. ‘pseudoharzianum' (Druzhinina 
et al. 2010b) what is also proposed for 
other pairs of sister species in the ge¬ 
nus Hypocrea/Trichoderma such as H. 
jecorina and T. parareesei (Druzhinina 
et al. 2010a, Atanasova, Jakutsch et al. 
2010) or T. longibrachiatum and H. ori- 
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entalis (Druzhinina et al. 2008). It is com¬ 
mon to argue against the assignment 
of spatial speciation modes to fungi 
as, although they may be distributed in 
overlapping ranges, they can occupy 
different microhabitats in those areas 
and therefore may be still spatially iso¬ 
lated (Burnett 2003). Our findings show 
the non-random presence of MOTUs 
attributed to both T. harzianum s. s. and 
H. ‘pseudoharzianum’ in at least four 
soil habitats. 


Biocontrol formulations 
may benefit from syn¬ 
ergistic action of highly 
opportunistic Hypocrea/ 
Trichoderma spp. 

A community of highly opportunis¬ 
tic Hypocrea/Trichoderma has been 
detected in previous cultivation-de¬ 
pendent studies: Migheli et al. (2009) 
showed that among 16 species iso¬ 
lated from highly disturbed non-rhizo- 
sphere soils in Sardinia (Italy) H. ‘pseu¬ 
doharzianum, T. spirale, T. gamsii. T. ha- 
matum or H. koningiopsis consistently 
co-occurred. A similar result was shown 
in the pioneering metagenomic study 
on Trichoderma in agricultural soils 
when representatives of Harzianum 
- Catoptron and Hamatum Clades, 
which cover species known for their 
antagonistic potential, were recov¬ 
ered from the same samples (Hagn et 
al. 2007). Zachow et al. (2009) applied 
metagenomic methods and tradition¬ 
al cultivation techniques to character¬ 
ize the diversity of fungi in rhizosphere 
of endemic plant species of Tenerife 
(Canary Islands) and showed the 
co-existence of extraordinarily high¬ 
ly antagonistic strains of Hypocrea/ 
Trichoderma. 

From a practical point of view, 
this demonstrates that the knowl¬ 
edge about infrageneric communities 
and interactions will be important for 
screening for Hypocrea/Trichoderma 
strains to be used for the biological 
control of soil borne plant pathogenic 
fungi. Several strains showing synergism 
with one another may be combined in 


certain biocontrol formulations, while 
on the other hand, the indigenous 
Hypocrea/Trichoderma should not 
have antagonistic properties against 
the introduced biocontrol strain(s). 

In this study we investigated one 
of the last remaining undisturbed eco¬ 
systems in Central Europe. However 
the diversity found in the national 
park largely resembles Hypocrea/ 
Trichoderma species composition 
in disturbed and agricultural soils 
(Migheli et al. 2009, Hagn et al. 2006). 
This finding indicates that the local 
highly opportuninistic species (from 
those listed above) should be among 
the major taxa screened for the best 
Trichoderma biocontrol strains. 


Both ITS rRNA and tefl 
phylogenetic markers 
have limited applicability 
for in situ diversity studies 
using the high-throughput 
methods 

Hagn et al. (2007) designed 
Hypocrea/Trichoderma-specific prim¬ 
ers for ITS1 fragment of the rRNA gene 
cluster. However, later studies showed 
that ITS1 is not sufficiently diagnostic as 
many species share the same allele. 
Meincke et al. (2009) set up genus-spe¬ 
cific primers with the reverse primer 
located in a still polymorphic and in- 
del-rich area of ITS2 30 bp upstream of 
the last genus-specific hallmark what 
makes several species undetectable. 
The six Trirev primers presented in this 
study amplify the entire diagnostic re¬ 
gion of ITS1 and 2 of all members of 
the genus. We have demonstrated 
the high specificity and selectivity of 
these primers as no MOTUs belonging 
to other fungi were recovered but the 
two novel alleles of ITS1 and 2 were 
detected. 

Hoyos-Carvajala et al. (2009) assert¬ 
ed the presence of paralogous cop¬ 
ies of ITS1 and 2 in some Hypocrea/ 
Trichoderma species. In order to test 
it we constructed a clone library for a 
randomly selected strain of T. asperel- 


lum (N=30) and found only one allele 
(data not shown) of ITS1 and 2 con¬ 
firming the absence of different alleles 
within a single genome. 

The applicability of ITS1 and 2 for 
larger metagenomic studies using 
high throughput sequencing methods 
remains questionable. On one hand, 
ITS-based quantification of species 
abundance in environmental samples 
will depend on the number of ITS cop¬ 
ies in the respective genome what 
may vary between different species. 
Furthermore, the growing number of 
Hypocrea/Trichoderma species share 
the same allele of ITS making these 
taxa indistinguishable by this locus 
(Samuels et al. 2006; Jaklitsch et al. 2006, 
Druzhinina et al. 2008; Atanasova, Jaklitsch 
et al. 2010). Another disadvantage of 
ITS1 and 2 comes from the fact that 
it is not appropriate for a design of 
species-specifc qPCR primers as its 
most diagnostic zones are surrounded 
by long mononucleotide stretches 
and/or have an intolerable GC con¬ 
tent (L. Bodrossy, I.S. Druzhinina, unpub¬ 
lished). Gazis et al. (2011) compared 
the three communities of endophytic 
fungi and revealed that ITS alone usu¬ 
ally underestimates the number of loci 
predicted by other loci. Sequences for 
the highly polymorphic 4 th large intron 
of tefl gene are also available for the 
majority of Hypocrea/Trichoderma 
species making it a good alternative 
to ITS1 and 2. This gene has a single 
copy in the genome and therefore 
becomes appropriate for quantita¬ 
tive assessments. Yet the current pub¬ 
lic database of Hypocreales tefl se¬ 
quences does not yet allow design 
genus-specific primers. In this study 
we attempted to develop species- 
specific tefl primers for those taxa, 
which we either detected among 
MOTUs or which we could expect to 
be present in temperate soils (based, 
for example, on Jaklitsch 2009, 2011). 
In vitro tests showed that only two out 
of 40 in silico designed oligonucle¬ 
otides are selective for target taxa 
while others demonstrated unspecific 
affinity to DNA of non-target species 
likely due to the polymorphic second¬ 
ary structures of the fragment. These 
results suggest that tefl intron alone is 
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also not appropriate for a large scale 
metagenomic analysis of the genus 
due to its hypervariablity. The applica¬ 
bility of more conserved markers such 
as rpb2 or chi 18-5 are currently tested 
in author's laboratory. 
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